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Abstract 

Photonic events with large missing energy have been observed in e + e _ collisions at a centre-of- 
mass energy of 189 GeV using the OPAL detector at LEP. Results are presented for event topologies 
consistent with a single photon or with an acoplanar photon pair. Cross-section measurements are 
performed within the kinematic acceptance of each selection, and the number of light neutrino 
species is measured. Cross-section results are compared with the expectations from the Standard 
Model process e + e~ — * vv + photon(s). No evidence is observed for new physics contributions to 
these final states. Upper limits on cr(e+e~ -> XY) ■ BR(X -> Y7) and cr(e+e~ -> XX) ■ BR 2 (X -> 
Y7) are derived for the case of stable and invisible Y. These limits apply to single and pair 
production of excited neutrinos (X = v* , Y = v), to neutralino production (X = X21 Y = x?) and to 
supersymmetric models in which X = Xx an( l Y = G is a light gravitino. The case of macroscopic 
decay lengths of particle X is considered for e + e~ — > XX, X — * Y7, when My « 0. The single- 
photon results are also used to place upper limits on superlight gravitino pair production as well 
as graviton-photon production in the context of theories with additional space dimensions. 
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1 Introduction 



We describe measurements and searches using a data sample of photonic events with large missing 
energy collected in 1998 with the OPAL detector at LEP. The events result from e + e~ collisions at a 
centre-of-mass energy of 188.6 GeV with an integrated luminosity of 177.3 pb -1 . The present paper 
builds on publications from earlier data samples at lower centre-of-mass energies |,§. This data-set 
at 189 GeV gives discovery potential in a new kinematic regime with about a four-fold increase in 
integrated luminosity. Measurements of photonic event production have also been made by the other 
LEP collaborations at centre-of-mass energies above the W pair threshold ||, including new results 
from L3 and DELPHI at yfe = 189 GeV @. 

The single-photon and acoplanar-photons search topologies presented here are designed to select 
events with one or more photons and significant missing transverse energy, indicating the presence of at 
least one neutrino- like invisible particle which interacts only weakly with matter. The event selections 
for these search topologies are similar to those used in our recent publication Jjj]. The single-photon 
search topology is sensitive to events in which there are one or two photons and missing energy which, 
within the Standard Model, are expected from the e + e~ — > 1/1/7(7) process^. The acoplanar-photons 
search topology is designed to select events with two or more photons and significant missing transverse 
energy which, within the Standard Model, are expected from the e + e~ — > 1/1/77(7) process. 

The single photon topology provides a direct measurement of the invisible width of the Z° and 
can probe charged and neutral triple gauge couplings. The acoplanar-photons topology can probe 
WW77 quartic couplings in the e + e~ — > i/ e zvy7 process. The neutral and quartic gauge coupling 
measurements will be described in forthcoming papers based on the event selections described herein. 

These photonic final-state topologies are sensitive to several different new physics scenarios. A 
generic classification is e + e _ — * XY or e + e~ — ► XX where X is neutral and can decay radiatively 
(X — * Y7) and Y is stable and only weakly interacting. For the general case of massive X and Y 
this includes conventional supersymmetric processes (X = X2,Y = Xi)- These topologies also have 
particularly good sensitivity for the special case of My ~ 0. This applies both to the production of 
excited neutrinos (X = z/*, Y = v) and to supersymmetric models in which the lightest supersymmetric 
particle (LSP) is a light gravitino and x± is the next-to-lightest supersymmetric particle (NLSP) which 
decays to a gravitino and a photon (X = x5,Y = G). The neutralino life-time in such models is a 
free parameter and so we also address the possibility of neutralino-pair production with macroscopic 
decay lengths. One type of new physics which could be seen in the single-photon topology is the 
production of an invisible particle in association with a photon. An example of this is graviton-photon 
production, e + e~ — * G7 This can occur within string theory models which allow gravitons 

to propagate in a higher-dimensional space but restrict Standard Model particles to the usual four 
space-time dimensions Q. Another type of new physics is the production of invisible particles tagged 
by initial-state radiation. One example^] is production of a pair of gravitinos, e + e~ — ► GG7, as in 
the superlight gravitino model ||. The acoplanar-photons search topology also has sensitivity to the 
production of two particles, one invisible, or with an invisible decay mode, and the other decaying 
into two photons. 

This paper will first describe the OPAL detector and the Monte Carlo samples used. A brief 
summary of the event selections will then be given, followed by cross-section measurements for e + e~ — * 
1/1/7(7) and e + e~ — * 1/1/77(7) an d comparisons with Standard Model expectations. The new physics 
search results will then be discussed. 

1 The photon in parentheses denotes that the presence of this photon is allowed but not required. 
2 The initial-state radiation diagram is one of many that contribute to this final state. 
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2 Detector and Monte Carlo Samples 



The OPAL detector, which is described in detail in Q, contains a silicon micro- vertex detector sur- 
rounded by a pressurized central tracking system operating inside a solenoid with a magnetic field of 
0.435 T. The region outside the solenoid (barrel) and the pressure bell (endcap) is instrumented with 
scintillation counters, presamplers and the lead-glass electromagnetic calorimeter (ECAL). The mag- 
net return yoke is instrumented for hadron calorimetry and is surrounded by external muon chambers. 
Electromagnetic calorimeters close to the beam axis measure luminosity and complete the acceptance. 

The measurements presented here are mainly based on the observation of clusters of energy de- 
posited in the lead-glass electromagnetic calorimeter. This consists of an array of 9,440 lead-glass 
blocks in the barrel (|cos#| < 0.82) with a quasi-pointing geometry and two dome-shaped endcap 
arrays, each of 1,132 lead-glass blocks, covering the polar angle^ range (0.81 < | cos^j < 0.984). Fully 
hermetic electromagnetic calorimeter coverage is achieved beyond the end of the ECAL down to small 
polar angles with the use of the the gamma-catcher calorimeter, the forward calorimeter (FD) and the 
silicon-tungsten calorimeter (SW). 

Scintillators in the barrel and endcap regions are used to reject backgrounds from cosmic ray 
interactions by providing time measurements for the large fraction (~ 80%) of photons which convert 
in the material in front of the ECAL. The barrel time-of-flight (TOF) scintillator bars are located 
outside the solenoid in front of the barrel ECAL and match its geometrical acceptance (| cos 6\ < 0.82). 
Tile endcap (TE) scintillator arrays [[!(]] are located at 0.81 < | cos 6 \ < 0.955 behind the pressure bell 
and in front of the endcap ECAL. 

The integrated luminosities of the data samples are determined to better than 1% from small- 
angle Bhabha scattering events in the SW calorimeter. Triggers [|i l] based on electromagnetic energy 
deposits in either the barrel or endcap electromagnetic calorimeters lead to full trigger efficiency for 
photonic events passing the event selection criteria described in the following section. 

The KORALZ fl2| and NUNUGPV98 fl3| Monte Carlo generators were used to simulate the 
expected Standard Model signal process, e + e~ — > vV + photon(s). For other expected Standard 
Model processes, a number of different generators were used: RADCOR [14 for e + e _ — > 77(7); 



BHWIDE p| for e+e" -» e+e"; TEEGG fjj for e+e" -> e+e'7; KORALW flj using grc4f |J 
matrix elements for e + e~ — > t~ vv{p[) and e + e _ — * vvqq, and KORALZ for e + e _ — > /x + /i~(7) and 
e + e~ — > t + t~(7). The Vermaseren program and grc4f were used for e + e~ — ► e + e~£ + £~. The 
expected contributions from each of these Standard Model processes were evaluated using a total 
equivalent integrated luminosity at least five times larger than the integrated luminosity of the data 
sample. 

To simulate possible new physics processes of the type e + e~ — > XY and e + e~ — * XX where 



X decays to Y7 and Y escapes detection, a modified version of the SUSYGEN |20] Monte Carlo 
generator was used to produce neutralino pair events of the type e + e~ — > X2X1 an d e + e~ — > X2X21 
X2 — ► X?7> w hh isotropic angular distributions for the production and decay of X2 an d including 
initial-state radiation. Monte Carlo events were generated at 48 (for XY production) and 42 (for XX 
production) points in the kinematically accessible region of the (Mx, My) plane. A Monte Carlo 



generator was written to simulate the superlight gravitino signature GG7, discussed in Section |5.1.4 . 



The KORALZ 1^7(7) sample was used to determine the efficiency for graviton-photon production in 
the context of additional space dimensions, by means of an event reweighting (see section 5.1.5| ). The 



3 In the OPAL coordinate system, 6 is the polar angle defined with respect to the electron beam direction and 
the azimuthal angle. 
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same procedure was also used to calculate the efficiency for GG7 production, and compared with the 
prediction of the direct Monte Carlo simulation. All the Monte Carlo samples described above were 
processed through the OPAL detector simulation [21|. 

Simulation of e + e~ — ► XX, X — > Y7 signal events with My ~ where Mx has a finite non- 
zero lifetime tx was implemented in the full simulation of the OPAL detector. In particular, the 
massive quasi-stable neutral particle X was assigned properties similar to a heavy neutrino for the 



purpose of propagation, and then propagated within the GEANT \22\ framework, according to its 
initial kinematics and tx before forcing the X — * Y7 decay. 



3 Photonic Event Selection 



This section summarizes the criteria for selecting single-photon and acoplanar-photons events. The 
kinematic acceptance of each selection is defined in terms of the photon energy, EL, and the photon 
polar angle, 9. In addition, the scaled energy, x 7 , is defined as Ey/Eb eam , and the scaled transverse 
energy, Xt, as x 7 sin#. 



Single-Photon - One or two photons accompanied by invisible particle(s): 

• At least one photon with x T > 0.05 and with 15° < 9 < 165° (| cos 6" | < 0.966). 
Acoplanar-Photons - Two or more photons accompanied by invisible particle(s): 

• At least two photons, each with x 7 > 0.05 and 15° < 9 < 165°, or one photon with EL > 
1.75 GeV and | cos 9 \ < 0.8 and a second photon with E 7 > 1.75 GeV and 15° < 9 < 165°. 

• The transverse momentum p7p of the two-photon system consisting of the two highest 
energy photons must satisfy p^V-^beam > 0.05. 



In each of the two cases we retain acceptance for events with additional photons in which the 
resulting photonic system is still consistent with the presence of significant missing energy. This reduces 
the sensitivity of each measurement to the modelling of higher-order contributions. Consequently, a 
large fraction of the kinematic acceptance of the acoplanar-photons selection is also contained in the 
kinematic acceptance of the single-photon selection. 



3.1 Single-Photon Event Selection 

The single-photon selection criteria are explained in detail in a previous publication [Q] . After defining 
the kinematic acceptance, additional cuts on cluster quality, forward energy, muon chamber and hadron 
calorimeter information, and multiphoton kinematics are used to remove cosmic ray backgrounds, 
beam related backgrounds and standard model physics backgrounds that have no missing energy. 
Then events are classed as either having or not having photon conversion candidates based on tracking 
chamber information. Each of these two classes has somewhat different additional selection criteria 
based on charged track activity, time-of-flight information, and other background suppression cuts. 

Several improvements have been made to the single-photon selection used in this analysis compared 
with the previous analysis [|lj . They are described in detail in Appendix A. 
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3.2 Acoplanar-Photons Event Selection 



The acoplanar-photons selection has classes of selection requirements similar to those of the single- 
photon selection. The details of the acoplanar-photons selection are described in our previous analysis 
[|]|; changes with respect to the previous analysis are minimal and are given in Appendix A. 

The requirement of a second photon in an event reduces many of the backgrounds which are 
otherwise a problem for the single-photon analysis. For this reason, the acoplanar photons selection 
can tolerate a lower energy threshold for the most energetic photon, as well as looser but more inclusive 
acceptance requirements for photon conversions. On the other hand, the single-photon selection has 
more acceptance for events having no time-of-flight information for the photons. In order to obtain 
the best overall acceptance for acoplanar-photons, we have added to the acoplanar-photons selection 
that part of the single-photon selection which contains two photons within the kinematic acceptance 
of the acoplanar-photons selection. This addition results in a relative increase in efficiency of 9.6% for 
Standard Model e + e~ —* 1/1/77(7) events. 



4 Selection Results 



The results of the single-photon and acoplanar-photons selections and the corresponding cross-section 
measurements and other measured event quantities are given below in sections [4.1| and 4.2. Results 
from both selections are summarized in Table [j]. 



4.1 Single-Photon 



After applying the single-photon selection criteria to the data sample, 643 events are selected. The 
expected contribution from cosmic ray and beam-related backgrounds is 4.6 ± 1.5 events. These back- 
grounds have been estimated from events having out-of-time TOF or TE information, but passing all 
other selection criteria, and from events selected with looser criteria that have been visually scanned. 
Of the expected physics backgrounds from plausible sources, e + e~ — > e + e~7, e + e~ — > £ + £~ vv(^\ 
e + e~ — ► e + e~£ + £~, e + e~ — > vvqq, e + e~ — ► n + fi~j and e + e~ — > r + r~7 have non-negligible contri- 
butionsQ. The total number of expected physics background events is 4.6 ± 0.5 and the contributing 
sources are summarized in Table |2[ The number of events expected from the Standard Model pro- 
cess e + e~ — > 1/177(7) as predicted by KORALZ is 679 ± 5 (stat) ± 14 (sys) where the systematic 
error is from experimental sources as discussed later. The number of events observed is consistent 
with the number expected from e + e~ — > 1/1/7(7) pl us the background. The efficiency for selecting 
e + e~ — ► 1/1/7(7) events within the kinematic acceptance of the single-photon selection is (82.1 ± 1.7)%. 
For both the single-photon and acoplanar-photons selections, efficiency losses due to vetoes on ran- 
dom detector occupancy range from about (3-5)%. Quoted efficiencies include these losses. The 
single-photon cross-section for ^/s = 189 GeV, accounting for detector and selection efficiencies and 
subtracting the estimated background, is 4.35 ± 0.17 (stat) ± 0.09 (sys) pb. 

The systematic error on the cross-section measurement from knowledge of the efficiency and nor- 
malisation is estimated to be 2.1%. The contributing uncertainties are summarized in Table 3. The 
dominant systematics are the selection efficiency uncertainty (1.5%) and the uncertainty on the de- 
tector occupancy estimate (1%). The event selection efficiency is controlled using a data sample of 

4 The expected number of events from the Standard Model process e + e~ — > vvvv^(^) is negligible and is neglected 
in both the efficiency and background estimates 
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around 1200 e + e~ — ► 77 events. In particular, the efficiency and time response of the TB and TE 
scintillators are measured from these data and small correction factors applied to the efficiency. The 
selection efficiency systematic includes the uncertainties on such corrections and an estimate of the 
residual uncertainty on the efficiency. The occupancy, estimated from random beam-crossing trig- 
gers, is typically 4% and we assign an error of 1%. We have compared the estimated efficiency using 
the NUNUGPV98 and KORALZ event generators. We find that the relative efficiency difference is 
(0.5 ± 0.4)% and assign an error of 0.5% to take into account the sensitivity to the modelling of 
the photon energy, angle and multiplicity distributions. This source of error is considerably reduced 
from our previous publications, largely because both generators that are now used are designed to 
model accurately 1/1/77 events. Additional systematic errors on the efficiency arise from uncertainties 
in modelling the material close to or inside the beam-pipe which accounts for photons which convert 
early (20% relative uncertainty) and from systematic effects in the track reconstruction and parameter 
estimation which affect the association of tracks to ECAL clusters. Uncertainties in the photon energy 
scale, resolution and angular measurement give small additional systematic contributions. 

The cross-section as a function of centre-of-mass energy is plotted in Figure [l]. Cross-section 
results from the current analysis and from our earlier publications [|], are plotted. The curve 
shows the predicted cross-section from the KORALZ event generator for the Standard Model process 
e + e~ — ► 1/7/7(7). The data are consistent with the prediction. 

Figure |2|a shows the recoil mass distribution, where the recoil mass M recoi i is defined as the mass 
recoiling against the photon (or against the two-photon system). The peak in the distribution at M% 
is due to a large contribution from the decay Z° — > uu. Figure ^b shows the polar angle distribution 
along with the e + e~ — > 1/1/7(7) Monte Carlo expectation. In both distributions, there is consistency 
between data and Monte Carlo, although we note a deficit in the radiative return peak and a slight 
excess in the low energy, high recoil mass region. 

The single-photon selection is designed to allow for the presence of a second photon in order 
to accept events from the e + e~ — * z/z/77 process. Thirty-six observed events are considered to be 
two-photon events (i.e. have a second photon with deposited energy exceeding 300 MeV and with 
15° < 6 < 165°), consistent with the expectation of 33.6 ± 1.5 events from the KORALZ Monte Carlo. 
Of these, 20 fall within the kinematic acceptance of the acoplanar photon selection, as compared to 
the KORALZ expectation of 23.5 ± 1.0. 

4.2 Acoplanar-Photons 

The acoplanar-photons selection applied to the data sample yields 24 events, in good agreement with 
the KORALZ prediction of 26.9 ±1.2 events for the Standard Model e + e~ — ► 1/1/77(7) contribution. 
The expected contribution from other Standard Model processes and from cosmic ray and beam 
related backgrounds is 0.11 ± 0.04 events. The selection efficiency for e + e~ — > 1/1/77(7) events within 
the kinematic acceptance of the selection is (66.4 ± 2.9)%. The corresponding cross-section is 0.204 ± 
0.043 pb, compared to a KORALZ prediction of 0.228 ± 0.002 pb. The OPAL measurements of the 
cross-sections at \/s = 130, 136, 161, 172, 183 and 189 GeV are summarized in Table [|. Results for 
y^s < 189 GeV have been taken from our previous publications fl],§]. 

The dominant source of systematic uncertainties is modelling of the reconstruction efficiency, espe- 
cially the simulation of the detector material and consequent photon conversion probabilities. Other 
sources arise from uncertainties on the electromagnetic calorimeter energy scale and resolution, on 
the integrated luminosity measurement, on detector occupancy estimates and from comparisons of 
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different Monte Carlo event generators for the process e + e — ► 1/1/77(7). The relative systematic error 
from all sources is 4.3%. 

The kinematic properties of the selected events are displayed in Figure [| where they are compared 
with the predicted distributions for e + e~ — > 1/1*77(7) obtained using the KORALZ generator nor- 
malized to the integrated luminosity of the data. Plot (a) shows the recoil mass distribution of the 
selected acoplanar-photon pairs (or of the two most-energetic photons in the case of events with three 
or more photons). The distribution is peaked near the mass of the Z° as is expected for contributions 
from e + e~ — > 1/1/77(7). The resolution of the recoil mass is typically 3-5 GeV for M reco ii ~ Mi. Plot 
(b) shows the distribution of the scaled energy of the second most energetic photon. Plot (c) shows 
the 77 invariant-mass distribution for which the mass resolution is typically 0.7-1.9 GeV. Plot (d) 
shows the distribution in scaled transverse momentum of the selected two-photon system. There is 1 
selected event having a third photon with deposited energy above 300 MeV and within the polar-angle 
acceptance of the selection. The corresponding expectation from KORALZ is 1.20 ± 0.08 events. 



5 Data Interpretation 

The results of the single-photon and acoplanar-photons selections are used to test the Standard Model 
and search for new physics contributions. 

For the XY and XX searches, we test the range of the following product branching ratios which are 
consistent with the data: cr(e + e~ -> XY) • BR(X -> Y7) and cr(e + e~ -> XX) • BR 2 (X -> Y7). Given 
lack of evidence for signal, we then set 95% CL upper limits on these quantities. This is done both 
for the general case of massive X and Y, and also separately for the special case of My ~ 0. All 
efficiencies are first evaluated under the assumption that the decay length of X is zero. For the XX 
search in the special case of My ~ 0, we evaluate the efficiency as a function of the non-zero lifetime of 
the X particle, thus quantifying the sensitivity of the search for the general case of non-prompt decay. 
For both the XY and XX searches, Monte Carlo samples were generated for a variety of mass points 
in the kinematically accessible region of the (Mx, My) plane. To set limits for arbitrary Mx and My, 
the efficiency over the entire (Mx, My) plane is parameterized using the efficiencies calculated at the 
generated mass points. As justified previously Q, we restrict the searches to Mx + My > Mi for the 
single-photon topology and to Mx values larger than about Mi/ 2 for the acoplanar-photons topology. 

5.1 Single-Photon 

With the single-photon topology we measure the number of light neutrino species and observe a rise in 
the cross-section at low photon energies which is consistent with the additional cross-section expected 
from charged current contributions to e + e~ — ► z/ e zvy(7). We also give results of searches for XY 
production, superlight gravitino pair production, and graviton-photon production in the context of 
models with extra space dimensions. 

5.1.1 Neutrino Counting 

Single-photon events are expected within the Standard Model from the process e + e~ — > 1/1/7(7). At 
the tree level, the cross-section for muon-type and tau-type neutrinos is attributable solely to s-channcl 
Z production with initial-state radiation. For electron-type neutrinos, the single photon cross-section 
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arises from Feynman diagrams corresponding to s-channel Z production with initial-state radiation 
and i-channel W exchange with radiation from the initial-state or the exchanged W. Higher order 
electroweak processes such as WW boxes are expected to give a negligible contribution to the cross- 
section. 

The results of the single-photon selection are used to measure the size of the s-channel Z production 
contributions and the W-related contributions (W amplitude squared plus the W-Z interference), 
which we parameterize in terms of Nj, and fw- The measurement of the pure s-channel Z production 
contributions is a direct measurement of the Z invisible width, which is related to the effective number 
of light neutrino generations, N„, defined as the ratio of the Z invisible width to the expected width 
in the Standard Model for one neutrino generation. The W-related contributions are parameterized 
by a multiplicative scale factor, fw, defined to be 1 for the Standard Model expectation. 

In order to measure N y and/or fw, we perform a binned fit to both the overall event rate and to 
the shape of the photon energy distribution by minimising the negative log- likelihood: 

- logL = - logVin^piE^^^, f w ) -> n obs (£ 7 )) (1) 

where n exp is given by 

n e xp(E-y, N„, fw) = fw n exp (E 7 ) + N u n exp (En f ). (2) 

Here n b s is the observed number of events, n exp is the expected number of events as a function of N v 
and fw derived by reweighting fully simulated 1/177(7) events, and V represents the Poisson probability 
for observing n„& s events given an expectation of n exp . n^ xp is the number of events for Standard Model 
W contributions^] and n^ xp is the number of s-channel events expected per neutrino generation; both 
nY xp and n^ xp were evaluated using the NUNUGPV98 generator^. 

The dominant systematic errors arise from the uncertainty on the selection efficiency and from the 
theoretical uncertainty on the expected number of events, with a minor contribution from the limited 
Monte Carlo statistics. The relative uncertainty on the selection efficiency is 2.1% and was discussed 
in Section ||. We assign a value of 2% to the theoretical uncertainty based on comparisons and 
estimated precisions of the NUNUGPV98, KORALZ and grci/^7 [^] event generators. In addition, 
an uncertainty related to the modelling of the photon energy spectrum is assigned, again estimated 
by comparing the NUNUGPV98 and KORALZ event generators. Other sources of systematic error, 
such as the uncertainty on the centre-of-mass energy and the expected background, result in negligible 
contributions. 

In a first step, we check the consistency with the Standard Model predictions by fitting for both Nj, 
and f w . The results are f w = 1.12 ± 0.13 (stat) ± 0.12 (sys) and N v = 2.63 ± 0.15 (stat) ± 0.11 (sys) 
with a correlation coefficient of —41%. Figure |] shows the 70%, 95% and 99% probability contours in 
the space of the two parameters, while Figure || shows the photon energy distribution for the events 
selected at 189 GeV, compared to the expectation with fw = 1.12, N v = 2.63. The data are seen 
to be in fair agreement with the Standard Model prediction for N v = 3 and fw = 1- In particular, 
the W contributions are observed with a high degree of significance, even when allowing N v to be 
unconstrained. 

In a second step, we assume the W contributions to be as predicted by the Standard Model and 
fit for N v . The result is: 

N v = 2.69 ± 0.13 (stat) ± 0.11 (sys). (3) 

5 Both the pure W t-channel, which is dominant, and its interference with the Z° s-channel are taken into account 
assuming Standard Model couplings of electrons and electron-type neutrinos. 

6 Using the NUNUGPV98 feature which allows i/giV/M an( A v y,^y.l{'y) events to be generated separately 
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Alternatively, one can assume N v = 3 as indicated by the precise, but less direct, measurements 
of the Z lineshape [24|, and fit for the relative size of the W-contributions. The result is fw = 
0.99 ± 0.11 (stat) ± 0.12 (sys), establishing that the W-contributions are observed and are consistent 
with the expectations from the Standard Model. 



5.1.2 Search for e+e - XY, X — »■ Y7 ; General case: My > 

We search for evidence of new physics processes of type e + e~ — > XY, X — > Y7. To select candidate new 
physics events, kinematic-consistency and degraded-resolution cuts are applied to events in the single- 
photon event sample. The kinematic-consistency cuts require that the energy of the most energetic 
photon be within the range kinematically consistent with mass values (Mx,My) after accounting 
for energy resolution effects. The degraded-resolution cuts reject events in which the most energetic 
photon is in the angular regions 0.78 < | cos 0\ < 0.82 or | cos 9\ > 0.95; energy resolution is significantly 
degraded in these regions and the Monte Carlo simulation is less reliable. These cuts are approximately 
94% efficient in selecting signal events within the kinematic acceptance of the single-photon selection, 
assuming uniform energy and cos# signal distributions. 

Two methods are used to search for, and place upper limits on, contributions from e + e _ — ► 
XY, X — ► Y7. The first is an event-counting method, based simply on the total number of events 
selected as candidates. The kinematic cuts, described above, are sufficiently loose that the results 
are relatively insensitive to the shapes of signal distributions in photon energy or cos 8 that result 
from the specifics of a particular model. Thus this method gives the more generally applicable cross- 
section limits for XY production. In many models, however, the distributions for the production 
and decay angles in the process e + e~ — > XY, X — ► Y7 are approximately isotropic. We therefore 
also perform a likelihood-based analysis to search for new physics contributions under the assumption 
of isotropically distributed production and decay angles. This results in a substantial increase in 
experimental sensitivity, since energy and angular distributions can now be used to differentiate signal 
from background. 

In both methods, Standard Model background is assumed to be from e + e~ — > 1/^7(7) only. Back- 
ground from other sources, including the estimated cosmic ray and beam-related background, is small 
and is neglected in the limit calculations. Since uncertainties due to energy scale and resolution effects 
in the low Mx — My region lead to large relative uncertainties in the estimated efficiencies, we restrict 
our search to the region in which Mx — My > 5 GeV. 



Event-Counting Method: 

Upper limits are calculated based on the total numbers of observed and expected background events 



using the method described in reference [25|. This procedure is similar to the method used in our pre- 
vious publication (T|, with the exception that only the kinematic-consistency and degraded-resolution 
cuts are applied in order to improve the generality of the results. 

The total signal efficiency within the single-photon kinematic acceptance is approximately 85%, 
varying by less than 1% over most of the (Mx,My) plane. However, the efficiency decreases sig- 
nificantly when the mass difference Mx — My becomes small, due to low photon energies resulting 
from the X — > Y7 decay. Contributions to the systematic error on the efficiency for selecting events 
from potential new physics sources are similar in nature and size to those discussed in our previous 
publication. The total relative systematic error varies from 2.5 to 5.5%, depending on Mx and My. 
These systematic errors have been treated according to the method in reference [pq] ; the effect on 
the upper limits is small. Uncertainties in the 1/1/7(7) background estimate have a more significant 
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effect on the upper limits, and are treated by means of a convolution within the limit calculations. A 
relative 3% uncertainty has been assigned to the expected background contribution, based on factors 
considered in the cross-section measurement and assigning a 2% theoretical uncertainty. 

There are 552 events selected after the degraded-resolution cuts are applied; the KORALZ expec- 
tation is 601 ± 14. The photon energy distribution for these events is shown in Figure ||. The number 
of selected events after further application of the kinematic-consistency cuts depends on the values for 
Mx & n d My; ranges for observed and expected events are displayed as contour plots in Figure @. The 
consistency between the number of observed events (n & s ) and the number of expected background 
events (rib) is given in Figure ||. Plotted is -Pfluct, the probability for observing at least n b s events 
given a background expectation of rib ± Ob events, defined by: 



-Pfluct — /I 



n=n oha 



7 



'b 



(4) 



This figure is intended to highlight regions in the (Mx, My) plane with an excess of events that could 
potentially be indicative of new physics. Such an excess would appear as a very low probability 
(-Pfluct ^ !%)• However, a deficit of events is actually observed for most of the (Mx,My) plane. 
This leads to high values for Pfl uc t; for example, for n b s 3> 1, a value Pfl uc t = 0.97 indicates a deficit 
consistent at approximately the 3% level. The calculated values of Pfluct range from 0.71% to 99.94%. 



Values with rather low and high probability occur partly because we consider several thousand 
points in the (Mx,My) plane; each of these has a different set of kinematic-consistency cuts which 
select different parts of the photon energy spectrum. The net effect is that almost every possible 
energy range is selected for some (Mx, My) value. For example, as seen in Figure ||, there is a deficit 
in the energy range from approximately 26 to 38 GeV; this range corresponds roughly to Mx = 113 
GeV and M Y = 74 GeV. 

The general structure of Figure || can also be understood from the photon energy spectrum. There 
is an overall deficit of approximately 49 events (~ 1.7<r). As seen in Figure ^, most of the deficit lies in 
the region of the radiative return peak, while there is an excess in the low-energy region. For regions 
in the (Mx, My) plane for which low-energy photons are not kinematically allowed, events populating 
the low-energy region of the energy spectrum are rejected and the significance of the remaining deficit 
is enhanced. In particular, this occurs when Mx — My is relatively large, and accounts for the large 
pfluct values in the low-My region of Figure ||. Conversely, regions with a small mass difference 
Mx — My correspond to low photon energies, which in turn correspond to the region of the photon 
energy spectrum in which there is an excess. This accounts for the small Pfl uc t values in the low 
Mx — My region of Figure || 

Cross-section times branching-ratio upper limits calculated using the event-counting method are 
shown in Figure || and range from 43 to 409 fb. Results for one Mx, My pair (Mx = 125, My = 63 
GeV) are shown as an addition to the expected 1/1/7(7) background in Figure [TJ. The added signal 
contribution is that which would be expected from a cross-section equal to the 95% CL upper limit. 



Likelihood Method: 

The likelihood-based analysis is a straightforward extension of an extended maximum likelihood 
fit 1 27]. Upper limits are calculated using information from the photon energy and angular distributions 
as well as the total number of observed events. The number of observed events, number of expected 
background events, and signal efficiencies are the same as in the event-counting method. 
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The likelihood function is given by 



L(a s ) = V{(n b + n s ) - n obs ) ■ Y[[f b P b (E t , 9$ + f s P s (E u 9 t )] 



(5) 



where the product is over all selected events in the data, and 

• n b = expected number of background events, 

• n s = e s a s C = expected number of signal events, with e s the efficiency for observing signal events, 
a s the signal cross-section times branching ratio, and C the integrated luminosity, 

• n obs = number of candidate events observed in the data, 

• V((n b + n s ) — ► n obs ) = e~^ rib+ns \n b + n s ) n ° bs / n obs \ = the Poisson probability to observe n obs 
events given an expectation of (n b + n s ), 

• P bjS (Ei, Oi) = probability density (normalized to one) for a photon i resulting from a background 
(b) or signal (s) process to have an energy Ei and polar angle 0j, and 

• h,s = n b,s/( n b + n s) = relative fraction of background or signal events. 

The expected number of background events, n b , is determined from a Monte Carlo sample of 
events generated by KORALZ. The signal efficiency, e s , is calculated by integrating an energy and 
angular distribution function over the region of kinematic acceptance (including kinematic-consistency 
and degraded-resolution cuts), and scaling the result to account for additional efficiency losses due 
to other cuts in the single-photon selection. This energy and angular distribution function serves 
as the the signal probability density function P S (E,9). The background probability density function 
P b (E,9) is obtained from a parameterization of the distributions formed by simulated vV^ip/) events 
generated by KORALZ and NUNUGPV98. In the region |cos#| < 0.72, the photon energy and 
angular distributions are independent, with the angular distribution given by 1/sin 2 9 and the energy 
distribution determined with a parameterization. In the region |cos#| > 0.72, energy resolution is 
dependent on cos 9. In this region, energy distributions are parameterized separately for slices of width 
0.1 in |cos0|; when properly normalized, these parameterizations together form the 2-dimensional 
probability distribution. 

The value for a s is not fixed; it is instead treated as a free parameter in the likelihood function. 
When properly normalized, the likelihood function can be thought of as a probability density function 
for a hypothesized a s consistent with the observed data. Restricting a s to be non-negative, the 95% 
CL upper limit a 95 is therefore determined from the following equation: 



The validity of the likelihood-based method was tested with Monte Carlo simulations. Two types 
of test were performed. The first test relies on the following definition of a 95% CL limit: if the true 
signal cross-section happened to be equal to the 95% CL limit, then at least 95% of a large number 
of identical experiments would result in data that are more signal-like (as defined by a likelihood 
comparison) than the experiment from which the 95% CL likelihood was derived. This type of test 
was performed for various combinations of large and small numbers of expected signal and background 



0.95 = 




(6) 
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events, using events generated randomly according to distributions derived from the expected XY 
signal and ^177(7) background. The second type of test involved the calculation of a large number of 
95% CL upper limits from samples generated with an identical signal cross-section. In this case, at 
least 95% of the calculated upper limits should be larger than the actual signal cross-section. Both 
types of test yielded agreement with expectations to within statistical errors (typically less than 0.5%). 
In addition, tests with fully simulated Monte Carlo events were used to test for any biases due to the 
signal and background parameterizations. No significant biases were found. 



Sources of systematic errors are the same as for the event-counting method, and are treated in the 
same manner. However, when calculating limits using the likelihood method, an additional 4% relative 
uncertainty in experimental sensitivity has been added to account for effects due to uncertainties in 
the background and signal parameterizations. The 4% estimate is the result of tests using different 
parameterizations; the uncertainty is again treated according to the method given in reference [26|. 



The 95% CL upper limits resulting from the likelihood method are shown as a function of Mx and 



My in Figure 10. The values range from 23 fb to 371 fb. 



5.1.3 Search for e+e - — > XY, X — > Y7 ; Special case: M Y « 



The case My ~ is applicable to excited neutrino models and to some supersymmetric models 
mentioned earlier. The results presented above include this case and no separate analysis is performed, 
although the results are highlighted here. The upper limits on <r(e + e - — > XY) • BR(X — ► Y7) for 
My ~ as a function of Mx range from 23 to 107 fb using the likelihood method, and from 43 to 
170 fb using the event-counting method. The limits are shown in Figure 11 for both methods. Also 
shown is the average limit expected in the absence of signal for the event-counting method. The actual 
limit is significantly lower than the expected limit (see discussion in Section 5.1.2j ). Table 5 gives more 
details of the search results. 



5.1.4 Gravitino Pair Production 



A supersymmetric model has been proposed in which the gravitino is very light ||. This model 
predicts a new source of single-photon events from the process e + e~ — > GG7. We use the single- 
photon topology to place constraints on the cross-section and therefore the gravitino mass in this 
model. The differential cross-section is given by 



d 2 a 



dxj d cos 9 



aGx* 
45 



m 



G 



(7) 



where a is the fine structure constant, Gn is the gravitational constant, mg is the gravitino mass, 



is the photon scaled energy (E^/E^ eSum ) and 9 is the polar angle, with 



f &6j (x,cos9) = 2(l-xf 



[1 - x)(2 - 2x + x 2 ) x(-6 + 6x + x 



leading to a soft photon energy spectrum. 
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(8) 



Two methods were used to determine the efficiency for observing events from e + e~ — * GG7. 
The first is by reweighting simulated 1/1/7(7) events generated by KORALZ. Each simulated event is 
weighted by /^(s, cos 9)/fKz(x, cos 9), where fxz is the 2-dimensional energy and angular distri- 
bution for the most energetic photon in events generated by KORALZ. We use a parameterization 
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Jkz which is valid up to photon energies of about 60 GeV, and thus we establish a kinematic accep- 
tance region for GG7 production: xt > 0.05, 15° < 9 < 165°, and E 1 < 60 GeV. We note that the 
cross-section for GG7 production with photon energies above 60 GeV is negligible. 

In addition, a Monte Carlo generator was written to generate events according to the distribution 
given in Equation ||. Initial state radiation was treated in a manner identical to that used by the 
EXOTIC p§ generator. Efficiencies were calculated and found to be identical to those using the 



reweighting method to within statistical errors. 

To achieve maximum sensitivity^, we place a maximum energy requirement on the observed photon 
energy: E~ < 30 GeV. With this requirement, we observe 195 candidates with an expected z^z/7 
background of 179.6 ± 5.4. The efficiency within the above kinematic acceptance region is estimated 
by the Monte Carlo to be 82.8%, and was assigned a relative 5% uncertainty. Using the event-counting 
method described earlier, we place a 95% CL cross-section upper limit of 293 fb, giving a lower limit^] 
on the gravitino mass of 8.7 fieV. 



5.1.5 Graviton-Photon: Search for Extra Space Dimensions 



There has been recent interest in string theory models which postulate the existence of additional com- 
pactified space dimensions; these models allow gravitons to propagate freely in the higher-dimensional 
space while restricting Standard Model particles to a 3+1 dimensional hypersurface 0. The funda- 
mental mass scale in this class of theories, Md, governs the rate of graviton production; it is possible 
that direct graviton-photon production could occur at significant rates at LEP2 energies with an ex- 
perimental signature of a single photon with missing energy (5|,^|. The single-photon search topology 
can therefore be used to place constraints on the fundamental mass scale Mo (or, equivalently, on the 
size of the extra dimensions) by placing limits on the graviton-photon cross-section. 

The differential cross-section is given by 

where 5 is the number of extra dimensions and Ss~i is the surface area of a <5-dimensional sphere of 
unit radius, with 
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x 4 cos 4 0]. (10) 



We use the event-counting method to place limits on graviton-photon production in the cases 
2 < 5 < 7. We use the same kinematic acceptance region as described in Section 5.1.4 : xt > 0.05, 

< 60 GeV, and 15° < 9 < 165°. The expected photon energy spectrum is soft, so in order to 
improve sensitivity, we require observed photon energies to be less than 34 GeV for all values of 8. 



The efficiency for observing events from e + e — > G7 is determined by reweighting simulated 1/^7(7) 



events generated by KORALZ in a manner equivalent to that described in Section 5.1.4: simulated 



7 The optimization condition chosen was that the expected upper limit on the cross-section for contributions from 
the e + e~ — > GG7 signal process be minimised, where the expected upper limit is defined as the average limit one would 
expect to set in the absence of signal. 

8 Evaluated with a — 



1 
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events are weighted by /g 7 (x, cos 0)//kz(x, cos 6), where fxz is the same KORALZ parameterization 
used earlier. The distributions fxz and fey have roughly similar shapes, and the calculated efficiency 
is insensitive to small differences in the KORALZ parameterization. 

We observe 208 candidates with an expected i/z/7 contribution of 196.0 ± 5.9 events. Assuming a 
5% relative uncertainty in the signal efficiency, upper limits on the cross-section and corresponding 
lower limits on Mjy for 2 < 5 < 7 have been calculated^. The results are given in Table 6 based on 
the convention^ of Equation 2 in reference || . 



5.2 Acoplanar-Photons 

5.2.1 Search for e+e~ -> XX, X -> Yj ; General case: My > 

The searches for e + e~ — > XX, X — ► Y7, both for the general case discussed here and the special case 
of My w discussed in 5.2.2, use the methods described in our previous publication Selected 
events are classified as consistent with a given value of Mx and My if the energy of each of the 
photons falls within the region kinematically accessible to photons from the process e + e~ — > XX, 
X — ► Y7, including resolution effects. The selection efficiencies at each generated grid point for 
the e + e~ — > XX, X — > Y7 Monte Carlo events are shown in Table [?| These values include the 
efficiency of the kinematic consistency requirement which is higher than 95% at each generated point 
in the (Mx, My) plane for which Mx — My > 5 GeV. Events from e + e~ — > 1/1/77(7) are typically 
characterized by a high-energy photon from the radiative return to the Z° and a second lower energy 
photon. The kinematic consistency requirement is such that the two photons must have energies within 
the same (kinematically accessible) region. As Mx and My increase, the allowed range of energy for 
the photons narrows, and fewer 2/1/77(7) events will be accepted. For the 24 selected events, the 
distribution of the number of events consistent with a given mass point (Mx,My) is consistent with 
the expectation from e + e~ — ► 2/2/77(7) Monte Carlo, over the full (Mx,My) plane. 

In our previous publication, because of concerns about the modelling of the Standard Model 
process e + e~ — ► 2/1/77(7), an limits derived from the acoplanar photons analysis were obtained without 
accounting for the expected background. The theoretical situation is now greatly improved, with two 
event generators [Til 13 agreeing to better than 1% for the total cross-section for this process within 



the kinematic acceptance of this analysis. For that reason, in this paper, all limits derived from this 
selection have been calculated taking the (KORALZ) background estimate into account. Figure ^ 
shows the 95% CL exclusion regions for cr(e+e- -> XX) • BR 2 (X -» Y7). The limits vary from 33 fb 
to 103 fb for M x > 45 GeV and M x - My > 5 GeV. In the region 2.5 GeV < M x - My < 5.0 GeV, 
the efficiency falls off rapidly (see Table |7|) . As this rapid fall increases the associated uncertainty in 
the efficiency, no limits have been set in this region. 

Systematic errors are due primarily to limited Monte Carlo statistics at the generated (Mx, My) 
points and the uncertainty on the efficiency parameterization across the (Mx, My) plane. The com- 
bined relative uncertainty on the efficiency varies from about (3-6)% across the plane (for Mx— My > 5 



GeV). All systematic uncertainties are accounted for in the manner advocated in reference 26]. This 
also applies to the limits for the My ~ case, presented in the next section. 

9 Evaluated also with a = j|g 
10 The convention of reference |i| differs. 
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5.2.2 Search for e+e - — > XX, X — »• Y7 ; Special case: My ~ 



For the special case of My ~ the kinematic consistency requirements differ from those used for 



the general case. One can calculate [29| the maximum mass, M^ ax , which is consistent with the 
measured three-momenta of the two photons, assuming a massless Y. A cut on M^ ax provides further 
suppression of the 1/1/77(7) background while retaining high efficiency for the signal hypothesis. This 
is discussed in more detail in reference We require that the maximum kinematically allowed mass 
be greater than Mx — 5 GeV, which retains (95.5j^' )% relative efficiency for signal at all values of Mx 
while suppressing much of the remaining 1/1/77(7) background. In our previous publication we also 
applied a recoil-mass cut at 80 GeV. In the case where background is not accounted for in the limit 
calculation, such a cut improves the expected sensitivity of the analysis. With background subtraction 
this is no longer the case. Therefore that cut has been removed. 



Figure 13 shows the expected M x iax distribution for signal Monte Carlo events with Mx = 90 
GeV and for e + e~ — * 1/1/77(7) Monte Carlo events. Also shown is the distribution of the selected data 
events. For the My ~ case, the efficiencies calculated from Monte Carlo events are shown in Table || 
after application of the event selection criteria and then after the cut on M™ ax . Also shown in Table |8] 
are the number of selected events consistent with each value of Mx as well as the expected number 
of events from e + e~ — > 1/1/77(7). The number of selected events consistent with a given value of Mx 
varies from 14, for Mx > 45 GeV, to 3 events at the kinematic limit. The expected number of events 
decreases from 15.8 ± 0.7 at M x > 45 GeV to 1.34 ± 0.07 consistent with M x > 94 GeV. 

Based on the efficiencies and the number of selected events, we calculate a 95% CL upper limit on 
o"(e + e~ — ► XX) • BR 2 (X — > Y7) for My ~ as a function of Mx- This is shown as the solid line in 
Figure [l4]. The limit is between 50 and 80 fb for Mx values from 45 GeV up to the kinematic limit. 
Also shown as a dashed line is the expected limit, defined as the average limit one would expect to 
set in the absence of signal. The limits can be used to set model-dependent limits on the mass of 
the lightest neutralino in supersymmetric models in which the NLSP is the lightest neutralino and 
the LSP is a light gravitino (X = xj,Y = G). Shown in Figure 14 as a dotted line is the (Born- 
level) cross-section prediction from a specific light gravitino LSP model [30] in which the neutralino 
composition is purely bino, with rrii R = 1.35m^o and mg L = 2.7m^o. Within the framework of this 
model, x? masses between 45 and 88.3 GeV are excluded at 95% CL. 

As described in section 2, the efficiencies over the full angular range have been calculated using 
isotropic angular distributions for production and decay of X. The validity of this model has been 
examined based on the angular distributions calculated for photino pair production in reference 



For models proposed in reference 32], the production angular distributions are more central and so 
this procedure is conservative. For a 1 + cos 2 9 production angular distribution expected for t-channel 



exchange of a very heavy particle according to reference [31], the relative efficiency reduction would 
be less than 2% at all points in the (Mx,My) plane. 



5.2.3 Search for e+e — > XX, X — > Y7 ; Special case: My ~ and macroscopic decay 
length 

As an extension to the special case of My ~ 0, we consider the sensitivity of the acoplanar-photons 
search to e + e~ — * XX, X — * Y7 when X has a macroscopic decay length. This extension evaluates the 
selection efficiency using signal Monte Carlo samples with various lifetime values tx- The X — > Y7 
decay is treated by a modified OPAL detector simulation package designed to handle delayed decays 
(as described in section ^) . The range of lifetimes considered extends from the near-zero lifetime value 
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of tx = 10~ 15 s to tx = 10~ 7 s at which point cr is 30 m. Such long lifetimes lead to decays outside 
the sensitive volume of the OPAL detector (determined by the outer radius of the ECAL), which 
implies a natural cutoff in sensitivity. 

As with the zero lifetime case, the (M^ ax > Mx — 5 GeV) requirement is applied. The relative 
efficiency of this cut decreases with increasing lifetime, but exceeds 90% for all (Mx, tx) combinations 
with high (greater than 50%) selection efficiencies. For (Mx,tx) combinations with lower selection 
efficiencies, the relative efficiency falls more rapidly with increasing tx, but does not drop below 
60% for the (Mx,tx) combinations considered. The loss in relative efficiency of the M£ iax cut is 
expected due to the definition of M^ ax , but the cut is maintained to allow comparison between the 
prompt decay and macroscopic decay length selection efficiencies. The resulting selection efficiencies 
for macroscopic decay lengths are given in Table ^. Also listed in Table ^ are the number of events 
observed in the data and the expected e + e~ — ► ^1777(7) background. The efficiencies indicate that 
sensitivity to macroscopic decay lengths is maintained up to lifetimes of 10 -9 s but falls rapidly for 
longer lifetimes. This is as expected, since for tx = 10 -9 s the decay length ranges from 30 to 3 cm 
for the values of Mx considered, so that the X — ► Y7 decays typically occur well within the confines 
of the detector. For lifetimes greater than 10 -9 s, the decrease in efficiency is two-fold. The primary 
loss is due to the increased decay length, implying fewer X — ► Y7 decays within the sensitive detector 
volume. Selection efficiency is also lost with increasing tx since larger lifetimes result in a delayed 
arrival time of the photon at the ECAL, which, if sufficiently late, causes the event to be vetoed due to 
the timing cuts imposed by the analysis. The effect of the event kinematics is also seen in the increase 
in selection efficiency for a given lifetime as Mx approaches its threshold value, corresponding to a 
drop in (3 from 0.84 at M x =50 GeV to 0.10 at M x =94 GeV. 

Contributions to the systematic error in the selection efficiency are (in order of significance) the 
limited Monte Carlo statistics, the efficiency parameterisation in the (Mx, Tx) plane and the modelling 
of the timing cuts. The total contributions are given in Table g 

Given the efficiencies and the number of selected events, the 95% CL upper limit on cx(e + e~ — * 
XX)-BR 2 (X -» Y7) for M Y w0 ciS Si function of Mx and tx is calculated, and is shown in figure |i~5|. The 
background-subtracted limits from the prompt decay analysis are used to define the profile as a function 
of Mx for tx = 10 -15 s. The interpolation of the limit to larger lifetimes is done by factoring the 
selection efficiencies of Tabl e ^ in to the 95% confidence level limits on a(e + e~ — * XX) • BR 2 (X — ► Y7) 



for My ~ 0. As in Section |5.2.2| the Born-level cross-section of the neutralino NLSP Gravitino LSP 
model of [30| is used to give an exclusion region in the (neutralino mass, lifetime) plane. This exclusion 
region is superimposed on Figure [H] and gives the excluded domain in the (Mx,tx) plane within the 
same specific light gravitino LSP model [pjj discussed above. 



6 Conclusions 



We have searched for photonic events with large missing energy in two topologies in data taken with 
the OPAL detector at LEP, at a centre-of-mass energy of 189 GeV. 

In the single-photon selection, which requires at least one photon with xt > 0.05 in the region 15° < 
6 < 165° (| cos 9 I < 0.966), 643 events are observed in the data. The background-subtracted cross- 
section measurement of 4.35 ± 0.17 (stat) ± 0.09 (sys) pb is consistent with the KORALZ prediction 
of 4.66 pb from the Standard Model 1/177(7) process. Interpreting the results as a measurement of 
the effective number of light neutrino species, we measure Nj, = 2.69 ± 0.13(stat) ± 0.11 (sys). We 
also observe significant W contributions to the cross-section with a rate consistent with the Standard 
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Model expectation. 

We calculate upper limits on the cross-section times branching ratio for the process e + e~ — * XY, 
X — * Y7 using two methods: an event-counting method, which is insensitive to energy or angular 
distribution shapes and is therefore relatively model-independent, and a likelihood method, which 
assumes isotropic production and decay angular distributions and has greater sensitivity. In the 
region of interest in the (Mx,My) plane, the limits vary from 43 to 409 fb using the event-counting 
method and from 23 to 371 fb using the likelihood method. These limits include the special case of 
My ~ 0, where the limit varies between 43 and 170 fb using the event-counting method and from 23 
to 107 fb using the likelihood method. We note that some of these limits are much more stringent 
than would be expected on average in the absence of signal contributions. 

We set a 95% CL cross-section upper limit on e + e~ — * GG7 production of 293 fb implying a 
lower limit on the gravitino mass of 8.7 /ieV in the super light gravitino model of reference pi. In 
the context of string theory models with extra space dimensions, upper limits on the cross-section for 
graviton-photon production of between 309 and 271 fb at 95% CL are set, giving lower limits on the 
fundamental mass scale varying between 1086 and 470 GeV for between 2 and 7 additional dimensions 
respectively. 

The acoplanar-photons selection requires at least two photons with scaled energy x 7 > 0.05 within 
the polar angle region 15° < 6 < 165° or at least two photons with energy > 1.75 GeV with 
one satisfying |cos#| < 0.8 and the other satisfying 15° < 9 < 165°. In each case, the requirement 
pJV-^beam > 0-05 is also applied. There are 24 events selected. The KORALZ prediction for the 
contribution from e + e~ — * 1/1/77(7) i s 26.9 ± 1.2 events; contribution from other sources is 0.11 events. 
The number of events observed in the data and their kinematic distributions are consistent with 
Standard Model expectations. We derive 95% CL upper limits on a(e + e~ XX) • BR 2 (X Y7) 
ranging from 33 to 103 fb for the general case of massive X and Y. For the special case of My ~ 0, 
the 95% CL upper limits on cr(e + e _ — > XX) • BR 2 (X — ► Y7) range from 50 to 80 fb. For the case of 
macroscopic decay lengths, these values range from 50 fb at short lifetimes to 8.8 pb for lifetimes as 
long as 10~ 7 s. 

The results of the acoplanar-photons search are used to place model-dependent lower limits on the 
Xi mass in a specific light gravitino LSP model |3(|. Masses between 45 and 88.3 GeV are excluded 
at 95% CL for promptly decaying neutralinos, while for a Xi lifetime of 10 -8 s (cr = 3m), masses 
between 45 and 81 GeV are excluded. 
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Appendix 



A Improvements to Event Selections 

Detailed descriptions of the event selection criteria are given in a previous publication This 
appendix describes changes and updates to the selections used in the present publication. 

A.l Single-Photon Selection 

There have been several improvements to the single-photon selection: 

• Non-conversion candidates in the endcap region (| cos#| > 0.82) are no longer required to have a 
good in-time associated TE hit, although candidates with associated out-of-time TE hits are still 
rejected. This change reduces sensitivity to the modelling of the material between the interaction 
point and the TE scintillators, and results in a relative increase in efficiency of about 19% in the 
endcap region with an associated increase in cosmic and beam-related background of about 1.3 
events. The special background vetos continue to be applied to all non-conversion candidates in 
the endcap regardless of the presence of TE timing. 

• The cluster extent cut has been modified for events in the endcap region. This cut had been fixed 
at lower centre-of-mass energies and was no longer fully efficient at = 189 GeV, particularly 
for beam energy photons. Events with the primary photon candidate in the region | cos 9 | > 0.82 
are now rejected if Acf) • sin 6 0/ In Ely > 0.05 radians^ where Acf) is the 4> extent of the cluster 
and Ey is the cluster energy in GeV. The result is a relative increase in efficiency of about 13% in 
the endcap region with no significant change in expected cosmic and beam-related background. 

• Small changes have been made to the timing requirements. The timing cuts for the photon 
candidates with TOF associated hits (|cos#| < 0.82) have been relaxed for certain periods of 
data-taking, resulting in an efficiency increase of 2.3% for events with a photon in | cos 6 \ < 0.82 
with a TOF associated hit. The expected cosmic ray background increases by about 10%. 

In addition, the majority of events rejected by poorly measured TOF times are recovered. In 
cases where the ECAL cluster and the TOF z-measurement from time difference differ by more 
than 40 cm, modified arrival times based on the ECAL cluster position and the time measurement 
at each end of the scintillator are constructed. Events are retained if either modified arrival time 
is within 5 ns of the expected time for a photon originating from the interaction point. As a 
result, the probability for a signal event to be rejected due to bad TOF timing falls from 0.6% 
to less than 0.1% with no measurable change in the cosmic and beam-related background. 

• The efficiency of the photon conversion consistency criteria has been improved substantially. In 
particular, all reconstructed tracks in the event (usually two for signal events) are now used to 
test association of a charged track to an ECAL cluster, and the association test also considers 
the disfavoured solution of the jet-chamber left-right ambiguity to account for cases where the 
incorrect solution has been chosen. 

The sin 6 dependence was parameterized for beam energy photons using samples of real and simulated e + e~ — > 77 
events. 
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• Improved redundancy in the rejection of beam-related backgrounds for conversion candidates 
with xt < 0.1 and | cos#| > 0.90 is introduced, based on the measured z-difference of the track's 
point of closest approach to the interaction point. 

• An additional beam-halo veto was implemented; events are rejected if ECAL clusters exceeding 
100 MeV in deposited energy are found within 20 cm in the r — (f> plane but on opposite ends 
of the ECAL endcap. This veto provides better rejection of beam-halo background and has 
negligible effect on signal efficiencies. 

• A minor addition was made to the special background vetoes used for non-conversion candidates. 
Events are now explicitly rejected if there is a large amount of hadronic energy in an HCAL 
cluster associated with the ECAL cluster of the primary photon candidate. This provides better 
redundancy with existing cuts in terms of background rejection, with negligible inefficiency even 
for beam energy photons. 

A. 2 Acoplanar-Photons Selection 

The acoplanar-photons selection is almost identical to that described in our previous publication. 
The only difference involves the timing requirements for candidates with associated TOF hits. These 
were relaxed for certain periods of the data-taking in the same manner as described above for the 
single-photon selection. 
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Quantity 


Topology 




Single-photon 


Acoplanar photons 


Nobs 


643 


24 


■^expected 


679 ± 5 ± 14 


26.9 ± 1.2 


N b kg 


9.2 ± 1.6 


0.11 ± 0.04 


Efficiency (%) 


82.1 ± 1.7 


66.4 ± 2.9 


O-meas(pb) 


4.35 ±0.17 ±0.09 


0.204 ± 0.043 


0"KORALz(pb) 


4.66 ± 0.03 (stat) 


0.228 ± 0.002 (stat) 



Table 1: Summary of results for the single and acoplanar photons selections. Shown are the num- 
ber of events observed in the data; the number of expected e + e _ — ► 1/1/7(7) (single-photon) or 
e + e~ — > 1/1/77(7) (acoplanar-photons) events based on the KORALZ event generator; the number 
of events expected from backgrounds; the selection efficiency within the kinematic acceptance of the 
selections (defined in section 3); the background-subtracted, measured cross-section within the kine- 
matic acceptance; and the expected cross-section based on the KORALZ generator. Errors when not 
marked are statistical and systematic added in quadrature. When two errors are shown, the first is 
statistical and the second is systematic. 



Background process 


Number of events 


e + e~ — > e + e~7 


0.3 ± 0.2 


e+e" -► t+tTvvij) 


1.80 ± 0.14 


e+e" -► e + e~£ + £- 


1.2 ± 0.4 


e + e~ — > visqq 


0.71 ± 0.09 


e + e~ — ► / u + / u _ (7) 


0.37 ± 0.08 


e + e~ — > r + r~(7) 


0.17 ± 0.05 


e + e~ — > 77(7) 


0.0 (< 0.05) 


Total physics background 


4.6 ± 0.5 


Cosmic background 


2.6 ± 0.6 


Beam related background 


2.0 ± 1.4 


Total background 


9.2 ± 1.6 



Table 2: Numbers of events expected from various background processes contributing to the single- 
photon event sample. The Standard Model background contributions are given by process. Also 
shown are the expected other backgrounds from cosmic rays and beam related sources. The errors 
shown are statistical. 
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0.7 


Trarlcinp' svstpmatirs 


0.5 


Event generator physics modelling 


0.5 


j;t scale and resolution 


0.4 


Integrated luminosity 


0.2 


Angular acceptance 


0.2 


MC statistics 


0.2 


Total 


2.1 



Table 3: Summary of the experimental systematic errors related to the efficiency and normalisation of the 
single photon cross-section measurement. 



^(GeV) 




«i 7 i!z(pb) 


130 


1.49 ±0.68 


0.626 ±0.010 


136 


1.23 ±0.56 


0.526 ±0.008 


161 


0.16 ±0.16 


0.330 ±0.018 


172 


0.32 ±0.23 


0.303 ±0.017 


183 


0.27 ±0.09 


0.247 ± 0.002 


189 


0.20 ± 0.04 


0.228 ±0.002 



Table 4: The measured cross-section for the process e + e~ — > 1/1/77(7), within the kinematic acceptance 
defined in section 3, for different centre-of-mass energies. For ^/s = 130 and 136 GeV the measurements 
are the weighted average of the results obtained from the 1997 data and the results obtained from 
the 1995 data. Results for y/s < 183 GeV are taken from our previous publications The final 

column shows the cross-section predictions from KORALZ. The quoted errors are statistical. 



M x 




E 

^max 


N obs 


■^expected 


-Pfluct 


iV 95 


(N 95 ) 


100 


25.1 


97.6 


390 


449 ± 13 


0.9921 


23.4 


52.6 


110 


30.5 


97.7 


371 


428 ± 13 


0.9917 


22.9 


51.0 


120 


36.3 


97.9 


358 


407 ± 12 


0.9832 


24.1 


49.5 


130 


42.6 


98.1 


342 


386 ± 12 


0.9768 


24.4 


48.0 


140 


49.3 


98.4 


320 


364 ± 11 


0.9812 


23.0 


46.4 


150 


56.5 


98.7 


293 


340 ± 10 


0.9899 


20.6 


44.6 


160 


64.0 


99.3 


258 


303 ±9 


0.9915 


19.0 


41.7 


170 


71.6 


100.1 


152 


182 ±5 


0.9844 


15.4 


31.5 


180 


78.8 


101.9 


15 


15.3 ±0.5 


0.5615 


9.3 


10.0 


185 


81.4 


104.0 


7 


6.2 ±0.2 


0.4213 


7.6 


7.1 



Table 5: Results for e + e — ► XY, X — > Y7, with My ~ 0. Shown are values for the mass Mx, the 
minimum and maximum energies (-^rnino Ernax) allowed by the kinematic consistency cuts, the number 
of observed events (N b s ), the number of events expected from e + e~ — > 1/1/7(7) (N ex pected)> the value 
of Pfiuct as defined in Equation ||, the calculated 95% CL upper limit on the number of new-physics 
signal events (iV 95 ), and the average value of iV 95 expected in the absence of signal ((iV 95 )). Upper 
limits are calculated using the event-counting method. Masses and energies are in GeV. 



24 



6 


e(%) 


& (fb) 


a 95 (fb) 


M D lower limit (GeV) 


2 


76.5 


241 


309 


1086 


3 


79.1 


233 


298 


862 


4 


81.4 


226 


290 


710 


5 


83.5 


220 


283 


605 


6 


85.4 


216 


276 


528 


7 


87.2 


211 


271 


470 



Table 6: Results for graviton-photon production in the context of extra dimensions. Shown are the 
number of extra dimensions 5, the efficiency e, the average upper limit on the signal cross-section 
(°"cxp) expected in the absence of signal, the measured 95% CL upper limit on signal cross-section 
cr 95 , and the corresponding lower limit on the mass scale Md- Efficiencies and cross-section limits are 



evaluated within the restricted kinematic acceptance described in Section 5.1.5 



M x 


M Y =0 


My = M x /2 


M Y = Mx - 10 


M Y = Mx - 5 


M Y = Mx - 2.5 


94 


70.3 ± 1.2 


70.7 ± 1.2 


64.1 ± 1.3 


38.9 ± 1.5 


5.1 ±0.7 


90 


69.4 ± 1.2 


71.9 ± 1.2 


60.4 ± 1.4 


40.5 ± 1.5 


4.6 ±0.7 


80 


71.5 ± 1.2 


69.3 ± 1.2 


59.7 ± 1.4 


39.5 ± 1.5 


4.7 ±0.7 


70 


71.0 ± 1.2 


70.8 ± 1.2 


60.4 ±1.4 


44.0 ± 1.5 


4.4 ±0.7 


60 


73.5 ± 1.1 


71.1 ± 1.2 


64.4 ± 1.3 


42.5 ± 1.5 


8.4 ±0.9 


50 


69.8 ± 1.2 


69.2 ± 1.2 


65.1 ± 1.3 


45.4 ± 1.5 


10.0 ±0.9 



Table 7: Acoplanar-photons selection efficiencies (%) for the process e + e~ — > XX, X — * Y7 at = 189 
GeV for various Mx and My (in GeV), after application of kinematic-consistency cuts. These efficiency 
values are used to perform the efficiency parameterization across the (Mx, My) plane. The errors 
shown are due to Monte Carlo statistics only. Efficiencies for the generated points at My = 20 and 
My = Mx — 15 are not shown, but are similar to those for My = and My = Mx/2. 



M x 


Selection efficiency 


Selection efficiency with 
Mg** > M x - 5 GeV 


Ndata 


■^V 1^77(7) 


94 


72.2 ± 1.2 


70.4 ± 1.2 


3 


1.34 ± 0.07 


90 


71.3 ± 1.2 


67.5 ± 1.3 


5 


2.40 ± 0.09 


80 


72.3 ± 1.2 


68.7 ± 1.4 


7 


4.81 ±0.13 


70 


71.4 ± 1.2 


69.2 ± 1.2 


9 


7.22 ±0.15 


60 


74.0 ± 1.1 


71.1 ± 1.2 


11 


10.05 ±0.18 


50 


70.2 ± 1.2 


67.7 ± 1.3 


14 


13.67 ±0.20 



Table 8: Acoplanar-photons event selection efficiencies (%), as a function of mass, for the process 
e + e~ — ► XX, X — > Y7, for My ~ at ^fs = 189 GeV. The first column shows the efficiency of the 
selection described in section 3.2. The second column shows the efficiency (%) after the additional cut 
on Mx ax - The third column shows the number of selected events consistent with the mass value Mx- 
The last column shows the expected number of events from the process e + e~ — > 1/^77(7) (KORALZ). 
The errors shown due to from Monte Carlo statistics only. 
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Selection efficiency with M™ ax > M x - 5 GeV 


lo gio( r x) 


M x =50 GeV 


M x =75 GeV 


M x =90 GeV 


M x =94 GeV 


-15 


65.2 ± 1.7 


68.2 ± 1.6 


71.1 ± 1.5 


72.1 ± 1.5 


-10 


63.3 ± 1.8 


67.8 ± 1.6 


72.2 ± 1.5 


71.4 ± 1.5 


-9 


59.9 ± 1.8 


71.7 ± 1.6 


71.3 ± 1.6 


72.1 ± 1.6 


-8.3 


25.8 ± 1.7 


43.5 ± 4.0 


49.7 ± 4.8 


51.4 ± 5.6 


-8 


13.9 ± 1.5 


25.8 ± 3.0 


33.0 ± 3.8 


37.8 ± 3.7 


-7.3 


2.5 ± 1.1 


2.7 ± 1.0 


5.6 ± 1.1 


10.7 ± 2.4 


-7 


0.8 ± 1.0 


1.0 ± 0.9 


2.2 ± 1.0 


5.1 ± 1.2 


Ndata 


14 


7 


5 


3 




13.67 ± 0.20 


5.88 ± 0.14 


2.40 ± 0.09 


1.34 ± 0.07 



Table 9: Acoplanar-photons event selection efficiencies (%), as a function of mass Mx, and lifetime 
tx for the process e + e~ — ► XX, X — > Y7, for My ss at = 189 GeV. The uncertainties in the 
selection efficiencies include the contributions due to the limited Monte Carlo signal sample size, the 
variation of the timing cuts, and the parameterisation of the fit efficiencies. Also included are the 
number of events found in the data and the number expected from the e + e~ — > 1/1/77(7) background 
that are consistent with the kinematic consistency cut for the Mx mass values considered. 
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Figure 1: The measured value of cr(e + e - — ► 7(7) + invisible particle(s)), within the kinematic 
acceptance of the single-photon selection, as a function of y/s. The data points with error bars are 
OPAL measurements at yfs = 130, 136, 161, 172, 183 and 189 GeV. The curve is the prediction for 
the Standard Model process e + e~ — * 1/177(7) from the KORALZ generator. 
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Figure 2: a) The recoil mass distribution for events passing the single-photon selection for the y/s = 
189 GeV data sample, b) The cos 8 distribution for the most energetic photon in the single-photon 
selection at ^/s = 189 GeV. In both plots, the points with error bars are the data and the histogram 
is the expectation from the KORALZ e + e~ — > vV^f^) Monte Carlo normalized to the integrated 
luminosity of the data. 
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Figure 3: Plots of kinematic quantities for the selected acoplanar-photons events for ^/s = 189 GeV. 
a) Recoil-mass distribution, b) Distribution of the scaled energy of the second photon (X2). c) 
Distribution of the invariant mass of the 77 system, d) Scaled transverse momentum distribution for 
the 77 system. The data points with error bars represent the selected OPAL data events. In each 
case the histogram shows the expected contribution from e + e~ — ► 1/1/77(7) events, from KORALZ, 
normalized to the integrated luminosity of the data. 
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Figure 4: The 70%, 95% and 99% confidence level contours in the fw^u plane resulting from 
the likelihood fit to the overall rate and the photon energy spectrum measured in the single-photon 
selection. The cross indicates the central values of the fit results, while the Standard Model expectation 
(N V =S, fw=l) is shown by the star. 
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Figure 5: Photon energy distributions for single-photon events. The points with error bars are the 
data. The solid histogram is the prediction for the values fw = 1-12, Nj, = 2.63 most consistent 
with the data. The dashed histogram is the expectation for the Standard Model values fw = 1> 
N„ = 3. The hatched region indicates the pure s-channel Z° contribution for N u = 2.63. All predicted 
distributions were calculated using the NUNUGPV98 generator. 
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Figure 6: The observed single-photon energy distribution after application of the degraded-resolution 
cuts. The data are shown as the points with error bars and are compared with the expectation from the 
Standard Model process e + e~ — ► vu r y('j), evaluated with KORALZ. The extra hatched contribution 
around 45 GeV illustrates the expected additional contribution from XY production which is excluded 
at 95% CL using the event-counting method for particular X and Y masses. The hatched contribution 
corresponds to Mx=125 GeV and My=63 GeV with a cross-section times branching ratio of 142 fb. 
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Figure 7: Contour plots for the y/s = 189 GeV data sample search for e + e~ — > XY, X — > Y7, 
after application of both degraded-resolution and kinematic-consistency cuts. For each set of mass 
values (Mx, My), a) shows the range of numbers of observed single-photon candidate events between 
each contour and b) shows the range of numbers of expected events from e + e~ — > 1/177(7) between 
each contour. The ^7(7) expectation was derived using KORALZ. Lines are drawn around the 
boundaries defined by M x + M Y = 189 GeV, M x = M Y , and M x + M Y = M z . The region in which 
Mx — My < 5 GeV is not considered. 
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Figure 8: Consistency between the number of observed events and the number of events expected 
from e + e~ — > 1/1/7(7) as estimated by KORALZ. Plotted is the probability for the number of expected 
events to fluctuate to the number of observed events or more, as described in the text. Values greater 
than 50% indicate a downward fluctuation where the number of observed events is fewer than expected. 
The boundaries and delineated regions are as defined for Figure [7|. 
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Figure 9: The 95% CL upper limit on cr(e + e" XY) -BR(X Y7) at y/s = 189 GeV as a function of 
Mx and My, using the event-counting method described in the text. The boundaries and delineated 
regions are as defined for Figure (?|. 
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Figure 10: The 95% CL upper limit on cr(e + e" -> XY)-BR(X -> Y7) at y/a = 189 GeV as a function of 
Mx and My, using the likelihood-based method described in the text. The boundaries and delineated 
regions are as defined for Figure 0. 
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Figure 11: The 95% CL upper limits on <r(e + e- -► XY) -BR(X -► Y7) at y/s = 189 GeV as a function 
of Mx, assuming My ~ 0, calculated using the likelihood-based method and the event-counting 
method as described in the text. Also shown is the average limit expected in the absence of signal, 
calculated using the event-counting method. 
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Figure 12: The shaded areas show 95% CL exclusion regions for cr(e+e- -»• XX) • BR 2 (X -► Y 7 ) at 
\fs = 189 GeV. No limit is set for mass-difference values Mx — My < 5 GeV, defined by the lower 
line above the shaded regions. The upper line is for Mx = My. 
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Figure 13: M^ ax distribution of the selected acoplanar-photons events (points with error bars). Shown 
as an unshaded histogram is the expected distribution from the Standard Model process e + e~ — * 
^77(7), evaluated using KORALZ and normalized to the integrated luminosity of the data. The 
shaded histogram shows the expected distribution for the signal process e + e~ — > XX, X — ► Y7 for 
Mx = 90 GeV with arbitrary production cross-section. 
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Figure 14: 95% CL upper limit on er(e + e~ -» XX) • BR 2 (X -» Y 7 ) for My ~ (solid line). Also 
shown is the expected limit (dashed line). The dotted line shows the cross-section prediction of a 
specific light gravitino LSP model |3(|. Within that model, Xi masses between 45 and 88.3 GeV are 
excluded at the 95% CL. These limits assume that particle X decays promptly. 
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Figure 15: The 95% CL upper limit on a(e + e~ -> XX) • BR 2 (X -> Y7) for My ~ as a function 
of Mx and log 10 (rx) (with tx in seconds), calculated using the efficiencies of Table |9| to extend 
the predictions of Figure Q to the case of X with a non-negligible decay length. Superimposed is 
the domain defined by the 95% CL exclusion limit (shaded region) under the assumption of a light 
gravitino LSP model [30]. 
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